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Using the general formalism for spherical gravitational collapse developed in 1 1], we investigate here the 
final fate of a spherical distribution of a matter cloud, where radial pressures vanish but tangential pressures are 
non-zero. Within this framework, firstly we examine the effect of introducing a generic small pressure in a well- 
known black hole formation process, which is that of an otherwise pressure-free dust cloud. The intriguing result 
we find is that a dust collapse that was going to a black hole final state could now go to a naked singularity final 
configuration, when arbitrarily small tangential pressures are introduced. The implications of such a scenario 
are discussed in some detail. Secondly, the approach here allows us to generalize the earlier results obtained on 
gravitational collapse with non-zero tangential pressure, in the presence of a non-zero cosmological constant. 
Finally, we discuss the genericity of black hole and naked singularity formation in collapse with non-zero 
tangential pressure. The treatment here gives a unified and complete picture on collapse final states, in terms 
of black hole and naked singularity formation, generalizing the earlier results obtained for this class of collapse 
models. Thus the role of tangential stresses towards determining collapse endstates emerges in a straightforward 
and transparent manner in our treatment. 
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I. INTRODUCTION 

Since the first proposal of the Cosmic Censorship Conjec- 
ture (CCC) [2], a great effort has been devoted in understand- 
ing the mechanism by which black holes and naked singulari- 
ties can form as the end state of the complete gravitational col- 
lapse of a massive object. However, no general proof, or even 
a suitable mathematically rigorous formulation of the CCC 
is available despite the efforts of many decades. Therefore it 
has become essential and inevitable that gravitational collapse 
within the framework of general relativity be studied in detail. 

The current situation is that, despite the great amount of 
efforts devoted in the past years in studying the CCC and its 
implications, the issue of final fate of complete gravitational 
collapse of a massive body is far from being entirely resolved. 
It is clear now how naked singularities have to be considered 
as a general feature of general relativistic physics, and that 
they may appear as the end-state of collapse in a broad variety 
of situations (see e.g. iH] and references therein). Still, it is 
not yet fully understood how different kinds of matter lead to 
the formation of either a black hole or a naked singularity. To 
give an example, how the presence of pressure in the form of 
tangential stresses and radial stresses might direct the collapse 
towards the formation of black holes or naked singularities is 
still a matter of debate. 

The first gravitational collapse scenario studied in detail 
within the framework of the general theory of relativity was 
the so called Oppenheimer-Snyder-Datt model @], that rep- 
resents a collapsing ball made of homogeneous pressureless 
dust. It is seen here that an event horizon forms before the 
appearance of the final spacetime singularity at the end of 
collapse, therefore hiding the singularity inside a black hole. 
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Further studies of inhomogeneous dust collapse though, such 
as those in the Lemaitre-Tolman-Bondi (LTB) models j5[], 
showed that naked singularities can indeed form during col- 
lapse, thus violating the cosmic censorship |6j]. However, 
the dust models are idealized scenarios that ignore pressures 
within the cloud. Therefore, collapse models with pressures, 
which are more physically realistic scenarios, were examined 
and it was soon discovered that for the collapse with non-zero 
pressures also, the end state is just as likely to be a naked sin- 
gularity, as for models without pressures (see e.g. |0]). 

In particular, collapse models with a vanishing radial pres- 
sure, but with non-vanishing tangential stresses, have been 
studied in some detail [8]. These models are of particular in- 
terest, because they provide a somewhat clear insight into the 
role of pressure towards determining the collapse final states. 

Using the formalism developed in [ 1 ] to study here the class 
of collapse models with a non-zero tangential pressures, three 
main results are obtained here. Firstly, we examine the effect 
of introducing small tangential pressure perturbations, on the 
black hole formation process, for a collapsing inhomogeneous 
dust cloud which is otherwise pressure-free. It is seen that 
a matter cloud that was collapsing to a black hole final state 
could change its course, and instead result in a naked singular- 
ity final configuration, when small tangential pressure pertur- 
bations are allowed in an otherwise pressure-free cloud. It is 
seen that the opposite scenario occurs as well. Secondly, we 
generalize here earlier results on gravitational collapse with 
pressure, in the presence of a non-zero cosmological constant. 
Finally, from the analysis of some collapse models, some in- 
teresting insights can be obtained on the genericity or other- 
wise of naked singularity formation in gravitational collapse, 
in the presence of a non-zero tangential pressure. Using the 
approach developed here, it is clearly seen how the black hole 
and naked singularity final states in collapse are distinguished, 
depending on the initial configurations of the matter density, 
the tangential stresses profiles and the velocities of the col- 
lapsing shells. Given this initial data, the Einstein equations 
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then fully determine the allowed evolutions that will lead to 
the formation either of a black hole or a naked singularity. 



II. EINSTEIN EQUATIONS AND REGULARITY 
CONDITIONS 



Various subcases of this class of collapse models have been 
studied in the past by different authors, as we point out below, 
but our treatment here gives a unified and complete picture 
on collapse final states in terms of black hole and naked sin- 
gularity formation, thus generalizing the earlier results. Thus 
the role of tangential stresses towards determining the final 
state of a collapsing matter cloud emerges transparently in our 
treatment. This helps us to see in a collective manner the dif- 
ferent results in this class of collapse models, which are oth- 
erwise somewhat scattered in a variety of approaches that do 
not always make it straightforward to see the overall conclu- 
sions in totality. Within this framework, the known results are 
easily recovered, and new scenarios and insights are obtained. 



In section [TT] we review as needed the general formalism 
for gravitational collapse in the presence of pressures, as re- 
stricted to the case of non-vanishing tangential pressure and 
vanishing radial pressure. The set of Einstein equations that 
governs the dynamical evolution of the collapsing cloud and 
the procedure to solve them is outlined. The regularity and 
physical reasonableness conditions that one would like to im- 
pose upon the matter cloud, in order to avoid any unphysical 
behaviour of the matter fields are described. The initial and 
boundary conditions that determine the dynamical evolution 
of the matter cloud are also discussed. Section [III] discusses 
the evolution of collapse, a function related to the tangent to 
the outgoing null geodesies at the singularity is evaluated, and 
it is the sign of this quantity that governs the final outcome 
of collapse in terms of either a black hole or naked singular- 
ity. The occurrence of trapped surfaces during collapse is dis- 
cussed, and the visibility of the singularity is studied through 
the analysis of outgoing radial null geodesies. 



The formalism is then applied in section[IV]to examine the 
effect of introducing small tangential pressures on the black 
hole formation process, for an otherwise pressure-free inho- 
mogeneous dust cloud collapse. In the section[V] we discuss 
the genericity of black holes and naked singularities obtained 
as collapse endstates, for this class of collapse models. The 
homogeneous and inhomogeneous dust cases are also briefly 
discussed, giving a review of known results for dust, obtained 
as a special case of the above treatment. In section IVII we 
generalize some earlier results on collapse with tangential 
pressure in the presence of a cosmological constant, which is 
asymptotically de-Sitter or anti-de Sitter spacetime. We study 
how the presence of a cosmological term in the equations af- 
fects the dynamical behaviour of the cloud. The study is then 
applied to the classes studied in section HVl by adding a cos- 
mological term and analyzing how the dynamics is altered. 
This constitutes a new class of solutions whose behaviour was 
not studied earlier. Finally, in section [VTIl we outline the in- 
teresting features of the above treatment in connection with 
previous results obtained by other authors. 



The general treatment to obtain the evolution and final fate 
of a spherical gravitational collapse with non-zero pressures 
was developed in and it constitutes the backbone upon 
which the present work is developed. One of the advantages 
of this approach is that it shows the final stages of collapse in 
a rigorous and yet transparent fashion, where all the crucial 
elements that enter towards the determination of the final out- 
come are enclosed in a single straightforward equation. It is 
clear that in order for the singularity to be visible, some out- 
going nonspacelike geodesies must come out of it to reach the 
faraway observers in the universe. This visibility of outgo- 
ing geodesies is shown to be related to the tangent to the null 
geodesies near the singularity, which in turn is shown to be de- 
pendent on the mass distribution, the pressure and the velocity 
profiles for the collapsing cloud. Once these quantities are set 
at the initial time, the evolution, as governed by the Einstein 
equations, determines entirely the final outcome of collapse. 

In spherical coordinates, the most general line element for 
a spherically symmetric collapsing cloud of matter depends 
upon three functions v, ip and R of the comoving radial coor- 
dinate r and the comoving time t, and it can be written as, 

ds 2 = _ e 2u(t,r) dt 2 + e 2mr) dr 2 + R ^ r) 2^2 _ (II. 1 ) 

We shall assume here through out that the radial pressure van- 
ishes. The energy-momentum tensor, expressed in co-moving 
coordinates then takes the form, 

It = -p; T; = 0; T e e = T* = Pe . (11.2) 

The Einstein equations then relate the three functions v,ip,R 
to the energy density and pressures of the system, and they 
can then be written as, 
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where the function F is the Misner-Sharp mass of the sys- 
tem and we have introduced the functions H(r, t) and G(r, t), 
which are defined as, 



H = e - 2v{r - v) R 2 , 

G = e ~mr,v) R >2 



(H.8) 
(H.9) 



Since there is a scaling gauge freedom available for the radial 
coordinate 7-, we can define without any loss of generality a 
scale so that at the initial time we have R(r, tj) = r. To this 
purpose, we introduce the scaling function v(r, t) such that, 
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The gravitational collapse condition is given here by R < 0. 
The area function R(r, t) can then be expressed as, 

R = rv (11.10) 

and the requirement of collapse is fulfilled once v < 0. 

The singularity is then achieved by reaching v = (which 
corresponds to R = 0), and it is realized along the singular- 
ity time curve t s (r). Note that as we approach the regular 
center of coordinates r = 0, when » ^ 0, the energy den- 
sity remains bounded, although the physical radius R goes to 
zero. Therefore, this indicates that the genuine curvature sin- 
gularity, where the physical radii of all the matter shells reach 
a vanishing value and the energy density diverges, is reached 
only in the case when v = 0. This is reflected in the scaling 
function by the requirement that v(r,t s (r)) = at the sin- 
gularity. To avoid the presence of shell-crossing singularities 
(that can be shown to be gravitationally weak and thus possi- 
bly removable |9| ). we further assume that, 

R'(r,t)>0. (11.11) 

The shell focusing singularity at v — is a true curvature 
singularity, and therefore is not removable from the spacetime, 
and its visibility or otherwise is the main issue that concerns 
us here. 

From equation (III. 3b , it follows that the Misner-Sharp 
mass F in this case does not depend on time (F = 
0), and therefore it must be a function of r only. We 
are then left with four equations in the six unknown 
p(r,t), PB {r,t), R(r,t), F(r), G(r,t), H(r,t). 

Collapse models with such a vanishing radial pressure, but 
non-zero tangential stresses were studied in the mass-area co- 
ordinates by Magli II 1011 . where the analytic structure of these 
solutions was derived. The algebraic equation which governs 
the nature of the central singularity was written and its depen- 
dence upon the initial configuration and the equation of state 
for the matter field was analyzed. Still, given the non-trivial 
nature of the problem, the task to study the outcome of col- 
lapse in full generality proved to be non-attainable. Further 
work by Goncalves, Jhingan and Magli [11] studied the final 
stages of collapse in the limited case of marginally bound col- 
lapse while the analysis of one of us [8] concentrated on mod- 
els with some restrictions on the pressure profiles (namely, 
assuming for the metric function v — v(R)). Also, the work 
by Harada, Nakao and Iguchi 111 211 analyzed the nakedness and 
strength of the central singularity obtaining some general con- 
ditions which were later specialized to the dust and 'Einstein 
cluster' scenarios. In this work the nakedness and curvature 
strength of the shell focusing singularity was studied in detail, 
providing a classification of the curvature strength of the sin- 
gularity based on some conditions (namely, the strong curva- 
ture condition, gravity dominance condition, and the limiting 
focusing condition). 

Another interesting model with vanishing radial stresses 
that has been studied is the so called 'Einstein cluster'. It 
is a specific subclass of models with tangential stresses that 
is constituted by a collapsing cloud of counter rotating parti- 
cles. Since the tangential stress is due to the particles' angular 



momentum, the 'Einstein cluster' is a particularly interesting 
model, because it reproduces the effects of rotation without 
departing from spherical symmetry. The analysis of the singu- 
larity forming in the 'Einstein cluster' IU3I1 shows once more 
how departure from dust models to include pressures does not 
forbid the formation of naked singularities. 

With the definitions ( IH.lOb . ( lH.8b , we substitute the un- 
knowns R, H with v, v. The definition of F(r) provides one 
of the two free functions and determines the energy density 
via equation ( III .41 . We can choose v as the second free func- 
tion and therefore evaluate pg through equation (III. 5b . which 
can be written in the form, 

Pe = \pR^ ■ (H.12) 

Without any loss of generality, we can obtain G from ( IH.61 l 
once we rewrite it as, 

G = 2v'^G (11.13) 
and we define a suitable function A(r, v) from 

A, v (r,v)=p'^ . (11.14) 
Then equation (III. 13b can be integrated to obtain, 

G{r,t) = b{r)e 2A{r ' v) . (11.15) 

The arbitrary function b(r) above that results from the integra- 
tion can be interpreted following the analogy with dust LTB 
models. Then b(r) turns out to be related to the velocity of 
the collapsing shells (more precisely to the velocity of the in- 
falling particles at the boundary of the cloud), and once we 
write it as 

b(r) = l + r 2 b (r) , (H.16) 

we can see that the values of bo = const, in the dust limit 
correspond to the bound (bo < 0), unbound [bo > 0) or 
marginally bound (bo = 0) Lemaitre-Tolman-Bondi models. 
Writing equation ( III. 71 ) in the form of a potential, 

fl 2 = e 2 " (I + G - l) =U(r,t) , (11.17) 

it is possible to analyze the dynamical behaviour of different 
pressure models with the usual phase space diagrams of clas- 
sical mechanics. Some examples will be discussed in section 
fVTl It turns out that depending on the velocity profile and the 
equation describing the tangential stresses, c ollap sing scenar- 
ios and bouncing scenarios are both allowed [ 10]. This marks 
a considerable difference with the LTB dust case where the 
potential U(r, t) is always decreasing, thus allowing only col- 
lapse to occur. 

Finally, from equation ill.li we get the equation of motion: 
e u [~F 

v = -V\Ir + g ~ 1 ' (IU8) 
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where the negative sign has been considered in order to de- 
scribe gravitational collapse. In some cases equation ( III. 18b 
can be integrated to obtain the function v(r, t), thus solving 
completely the system of Einstein equations. In other situa- 
tions, considerations on the behaviour of v(r, t) near the cen- 
ter of the collapsing cloud suffice to provide the analysis of 
the dynamical behaviour of the system approaching the sin- 
gularity. 

A viable physical description of the gravitational collapse 
of a massive body must be subdued to certain constraints and 
regularity conditions in order to ensure that the physical quan- 
tities that appear in the model are well defined and regular 
everywhere. 

Requiring regularity of the energy density at the initial sur- 
face from which the collapse commences, and at the center of 
the cloud throughout collapse, imposes some constraints on 
the classes of mass profiles that can be considered. Specifi- 
cally, the mass profile F(r) must have the form, 

F(r) = r 3 M(r) , (11.19) 

which immediately implies the regularity of the density func- 
tion at the center of the cloud at all regular epochs. By this 
prescription, we ensure that the Misner-Sharp mass F(r) van- 
ishes at the center of the cloud and that the energy density p 
away from the singularity has a finite value at the center. A 
further condition that may be imposed sometimes is that the 
energy density does not present a cusp at r = 0, and that it 
is a smooth function of the radial coordinate r. This can be 
achieved by requiring that the linear term of the energy den- 
sity expansion near the center vanishes, which is reflected in 
the requirement that M'(r) = at r = 0. Of course, this 
is an additional requirement and in a general physically real- 
istic scenario, this may or may not be realized. In any case, 
requiring smoothness brings in some additional mathematical 
simplification. 

The requirement that there is no pressure gradient at r = 
implies that the tangential pressure must become equal to 
the radial pressure (i.e. zero in our case) at the center. This 
imposes the condition that pg goes as r at least, near the center, 
and implies regularity for the metric function v via equation 
(1H.51 >. This condition then implies that we have in general, 

v(r,v) = r 2 g(r,v) , (11.20) 

with g being a regular function and u(r,t) = v(r,v{r,t)). 
Regularity of v implies regularity of the function A via equa- 
tion (1H.14I I. and we have, 

A(r,v) =r 2 a(r : v) , (11.21) 

with a being again a regular function. The relation between a 
and g is given by the equation (1H.14I I. which can be written as 

2ff + rg' 

a, v = — -jjj — • (11.22) 

As we have seen before, a specific behavior of R is required 
in order to avoid shell crossing singularities. The choice of 



R satisfying R' > implies further regularity for the en- 
ergy density. This in turn gives some constraints on the non- 
vanishing terms in the expansion of the energy density at a 
constant time surface. Namely, we find that in order to avoid 
the shell crossing singularities, we must require M' < 0, 
which, from a physical point of view, is a reasonable con- 
dition since we expect the density to be higher at the center of 
the cloud, and decreasing as we move away from the center. 
In the case when M'(0) = 0, p turns out to be a decreasing 
function from the center outwards if M"(0) < 0. 

As noted, from the regularity of v at the center of the cloud, 
we get the behaviour of b as in equation dll.161 ). and this turns 
out to be related to the velocity profile, and thus the kinetic en- 
ergy, of the infalling particles. Positivity of G is thus reflected 
in the condition that 

b a > -i . (11.23) 

Furthermore, the requirement that R' > implies that moving 
away from the center, the contribution of the kinetic energy to 
the total energy of the collapsing cloud decreases |9|. There- 
fore, the requirement of avoidance of shell crossing implies 
that some further condition must be imposed upon bo(r). In 
general, from the Misner-Sharp mass we get 

F'{l-G+H) >F{H'-G') , (11.24) 

which, in the limit of the dust case, becomes, 

b M' > Mb' . (11.25) 

The marginally bound case (bo = 0) is obviously allowed, 
and we see that considering only the constant values for the 
function ba (r) (which correspond in the dust case to the bound 
and unbound LTB models), the bound case, i.e. bo < might 
imply the occurrence of shell crossing if M' > 0. The bound 
case is therefore allowed under the restrictions imposed as 
above. We note from equation ( 1H.231 I that the permitted val- 
ues for bo have a limited range, namely bo € ( — ^, 0), where 

r b 

rfc is the boundary of the cloud (an infinite cloud cannot there- 
fore have a bound velocity profile). On the other hand, the 
unbound models, which are also allowed, describe a collapse 
scenario that is not uniquely gravitational. Some additional 
energy is impressed to the infalling particles from outside, 
causing the particles to have a non zero velocity at all times. 
We note that in the unbound case, solutions with R = at the 
initial time are not possible. 

In the more general situation where bo — bo (r), the fact that 
certain velocity profiles are allowed or not might depend on 
the boundary conditions imposed on the matter cloud. There- 
fore, for some specific cases, both positive and negative func- 
tions bo (r) can be considered, without causing the presence of 
shell crossings. 

At this point, we note that the requirement that the regular- 
ity conditions must be satisfied constrains the class of allowed 
models and the allowed equation of state. Regularity would 
be reasonable to demand from the point of view of the phys- 
ical validity, since it guarantees that the metric functions and 
physical quantities are well behaved at any regular epochs of 
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collapse, and at the center. We see that this imposes certain 
non-trivial restrictions on the possible choice of the matter 
models. In fact, we notice how regularity and the choice of 
a linear equation of state of the type pg = kp, as was assumed 
in some previous works (see for example, 111410 . are mutu- 
ally exclusive. This is immediately seen from equation ( III. 5b 
which now becomes, 

/?' 

v' = 2k—, (11.26) 
R 

and once integrated and evaluated at the initial time t = t j it 

gives, 

Vi{r) = 2fclnr + const. (11.27) 

Thus the metric function diverges as r tends to zero, while on 
the other hand regularity demands that v ~ r 2 at all times, 
near the center. We thus see that in comoving coordinates the 
assumption of a linear equation of state leads to a breakdown 
of the coordinate system at the center, or the metric functions 
not being regular. This would imply that the initial data set 
is not regular, or that a more suitable coordinate system is 
required. The inspection of the Kretschmann scalar K for this 
model shows that for any k ^ 0, the central shell presents a 
curvature singularity (where K diverges) at all times, and thus 
the initial data surface is not regular. This is confirmed also by 
the divergence at r = 0, of the energy density p, as evaluated 
from equation ( IH.41 i. It is seen that at the initial epoch also, 
the energy density diverges necessarily at the center. For the 
tangential pressure models, the choice pe — kp turns out to be 
consistent with the regularity condition F(r) = r 3 M (r) only 
in the dust case, where k = 0. 

In order for the matter model to be physically realistic, suit- 
able energy conditions, in particular the weak energy condi- 
tions, must be satisfied. That is, the energy density as mea- 
sured by local observers must be non-negative, and for any 
timelike vector V a we require, 

T ab V a V b > , (11.28) 

which reduces to 

P > 0, p + pe > . (11.29) 

Since the physical radius R is always positive, to ensure the 
positivity of energy density we must require F' and R 1 or their 
ratio to be positive. As we have seen, positivity of R' ensures 
that there are no shell crossing singularities occurring, while 
the fact that F' is positive is granted by the choice of positive 
M{r) with a suitable behaviour. From equation ( IH.12I I we see 
that positivity of p + pg is then given by 

\p > -i <n.30) 

and is satisfied automatically (and not only) by any increasing 
function v. The sign of v 1 is the only factor deciding positivity 
or negativity of the tangential pressure and it is worth noting 
that negative tangential stress profiles are also allowed and 



they satisfy energy conditions provided that the relation ex- 
pressed in equation (III. 3 Oi l is satisfied. Furthermore, for small 
values of r, from equation (III.301 > we see that regardless of 
the values taken by v, there will always be a neighborhood 
of r = for which \pg\ < p and therefore p + pg > is 
automatically satisfied in a close neighborhood of the center. 

Finally, specifying the initial conditions for the collapsing 
matter cloud consists in prescribing the values of the three 
metric functions and that of the density and stress profiles at 
the initial time U, given as functions of r lfl5ll . These can be 
written as, 

p{r,U) = pi(r), pg(r,U) =p 6i (r), 

R{r,U) = Ri(r), u(r,v(r,ti)) = Vi(r), ip(r,U) = ipi(r). 

As noted earlier, at the initial time we have taken the scale 
such that Ri = r, which implies that the scaling function 
v = 1 at the initial epoch. Furthermore, from R^ = 1 we 
get v'(r,ti) = 0. 

Not all of these initial value functions can be chosen arbi- 
trarily. In fact, the scaling function, the choice of the mass 
profile, and the Einstein equations impose some constraints. 
From equations ( III. 4b and (III. 19b , we see that 

Pi = = 3M (r) + rM'(r) , (11.31) 

while from equation (III. 12b we get, 

Pe, = \v'i{r)Pi{r)r . (11.32) 
Also, from equation ( III. 201 ) we can write, 

Vi{r) = r 2 g{r, v(r, t 4 )) = r 2 9l {r) . (11.33) 

In turn, v± can be related to the function a, defined by equation 
( III. 211 ). at the initial time via equation dH.14b . 

a. v (r, v(r, U)) = (a. v ) l = 2g t + rg[ . (11.34) 

Finally, we can relate the initial condition for ip to the initial 
value of the function A(r, t) through equation ( III. 9b and equa- 
tion dll.141 ). 

A{r,v{r,U)) = Mr) = - ^lnb . (11.35) 

If we want the initial collapse configuration to be not 
trapped, in order to see how different initial matter configu- 
rations evolve forming trapped surfaces, we must further re- 
quire, 

^ = r 2 M(r b ) < 1 , (11.36) 

which will impose some restrictions on the possible choices of 
the radial boundary in accordance with M(r). This condition 
reflects on the initial configuration for G and H since ^ < 1 
implies H < G. Therefore to avoid trapped surfaces at the 
initial time the velocity of the infalling shells must satisfy 

— R> Vbe A+v . (11.37) 
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This means that the scenario where collapse commences from 
an equilibrium configuration where R = 0, can be taken only 
as a limiting case. 

As we noted above, considering only the presence of the 
tangential pressures implies conservation of the Misner-Sharp 
mass during the collapse. Therefore, matching to an exte- 
rior spherically symmetric vacuum solution leads inevitably to 
consider the Schwarzschild metric [16]. The collapsing cloud 
has a compact support within the boundary taken at r = r^, 
and the pressure of the matter composing the collapsing star 
is assumed to vanish at the boundary. Matching conditions 
imply continuity of the metric and its first derivatives across 
this surface, although it should be noted that eventual discon- 
tinuities in the first derivatives of the metric across the bound- 
ary may be interpreted as an additional distribution of matter 
acting as a layer separating the interior from the exterior. In 
this case, the jump of the second fundamental form across the 
boundary is related to an additional stress-energy tensor layer 
confined to the boundary. Such a matching is in principle al- 
ways possible, and together with the regularity conditions and 
the energy conditions, might impose some further restrictions 
on the initial configurations allowed. 



III. COLLAPSE EVOLUTION AND FINAL STATE 

As we show below, the possible future evolution of the col- 
lapsing matter cloud depends upon the choice made for the 
matter content, namely the two free functions that determine 
the mass profile and the tangential pressure profile, and the 
initial conditions. 

In order to investigate the final outcome of collapse, we will 
study the time curves that lead the shell of matter labeled with 
a specific coordinate value r to fall into the singularity. As 
already mentioned, the spacetime singularity is achieved for 
v = 0, and the divergence of the energy density (and eventu- 
ally of the tangential pressure) indicates that it must be a true 
curvature singularity. 

Given the fact that the function v(r, t) is monotonically de- 
creasing in time, we can invert it and consider the time t as a 
function of the variables r and v. This is equivalent to consid- 
ering the area-radius coordinates (r, R) as done in some ear- 
lier works (see e.g. lflol[l2ll ') with the relation between R and 
v given by equation dll.lOb . The function t(r, v) will identify 
the time at which the shell labeled by a specific value r reaches 
the event v. In this manner, the occurrence of the singularity 
is described by the time curve t s (r) — t(r,0). The tangent 
of the singularity curve t s (r) at r = 0, i.e. at the central sin- 
gularity, is then seen to be related to the eventual existence of 
outgoing radial null geodesies escaping away from the singu- 
larity, and therefore to the possible visibility or otherwise of 
the singularity itself to faraway observers. 

From equation ( III. 181 we can write, 



dt = - 



I MM , b(r)e 2A (r,v) _i 

v ~~ ~ H p — 



dv 



(III.l) 



to describe the collapse. Integrating the above equation, while 
treating the radial coordinate rasa constant, we get 



t(r,v) = U + 



M , be 2A -l 
— T" 2 



-.dv 



(III.2) 



It follows therefore that, the time taken by the shell at any 
given value of r to reach the singularity, located at v = 0, is 
given by, 



t.(r) = U + 



M + be 2A -l 



-.dv 



(III.3) 



We note that t(r, v) is in general a differentiable function at 
least outside the singularity, which would be the case because 
of the regularity of all the functions involved. That is the case 
in many of the known collapse models. In any case, when- 
ever it is at least a C 2 function, we can always write it as an 
expansion near the central shell r — in the form, 



t(r,v) =t(0,v) +r X i(v) +o(r 2 ) , (III.4) 



with 



t(0, v) = U + 



-.dv (111.5) 



v M(0) 



■&o(0) + 2a(0,i>; 



and Xi(v)= | r=0 , given by, 



Xi(v) 



1 M' , rb'e A +rbA'e A -2(be A -l) 



( M 

V v 



be 2A - 
— 73- 



dv . 
= ° (111.6) 

Note that the above quantities are defined only if A1 ^ ' + 
&o(0) + 2a(0,v) > 0. This may be called the reality con- 
dition for the gravitational collapse to take place. The choices 
of b(r) and A(r, v), made from dII.16l > and ( 111.21 1 ) imply that 

as r — > 0, we must have he r ^ — > bo + 2a(0, v), and therefore 
we see how v cannot have constant terms or terms linear in r 
in a close neighborhood of r = 0. This justifies the regularity 
conditions for b(r) and v(r) that were deducted earlier. Fur- 
thermore, we can see how the regularity condition for A, that 
comes directly from that for v, ensures that t(0, v) and \i do 
not diverge. Evaluating xi explicitly we get 

m 1 f 1 zgjtW + MM , mT7 , 

Xi( u ) = -7t/ - — —dv. (IH.7) 



5 (^+b (0) + 2a(0,v))' 



with the minus sign chosen in accordance with v < 0, in order 



By continuity, we can therefore take that the singularity 
curve f s (r), which corresponds to v = 0, is differentiable, 
or in any case even if it is only a C 2 function, we can take 
it that it has at least a well-defined tangent at the singular- 
ity. This is the case in important known collapse models such 
as the Tolman-Bondi-Lemaitre dust collapse, or the 'Einstein 
cluster' models, and also for the Oppenheimer-Snyder-Datt 



7 



homogeneous dust collapse. The singularity curve is com- 
puted as the time it takes for the shell located at r to reach the 
singularity and can then be written as, 



t s (r) = t Q + r X i(Q) + o(r 2 ) 



(III.8) 



where to = t s (0) is the time at which the central shell be- 
comes singular. It can be shown that for any generic initial 
configuration there exist classes of dynamical collapse evolu- 
tions that have a well defined singularity curve. 

While the collapse evolves, the trapped surfaces may form 
in the collapsing cloud. If the trapped surfaces form at a time 
anteceding the formation of singularity then the latter will be 
covered, thus giving rise to a black hole at the end of the col- 
lapse. If, on the other hand, the trapped surfaces do not form, 
or form at a later time after the singularity, then the singularity 
might be visible to faraway observer 01711 . 

The boundary of the trapped region is delimited by the ap- 
parent horizon. The equation for the apparent horizon in a 
spherically symmetric spacetime is given by g %] R^Rj = 0, 
which for the metric (III, lb becomes G — H = 0. From equa- 
tion (III. 7b we can write it as, 



F 
R 



(111.9) 



The behavior of the apparent horizon for the collapsing cloud 
with vanishing radial pressure is analogous to that of the dust 

case, since ■£ = r M ^ in both these cases. If we note that 
as v — > we must have r — > on the apparent horizon, it 
is immediately clear that the only singularity that can eventu- 
ally be visible in the present case is that at the center of the 
collapsing cloud. In fact, since it is not possible to satisfy 
1 — ^ > near the singularity but away from the center, 
the region surrounding the singularity cannot be timelike, and 
therefore any singularity that might eventually form near the 
center with r > must not be visible. 

The difference with the dust case is then given by the dif- 
ferent possible behaviours of the function v(r,t). It is this 
function, as given from equation ( III. 18b . that determines the 
apparent horizon curve in dependence of the mass, velocity 
and the tangential stress profiles. We can see the apparent 
horizon as a curve r a h (t) given by 



Vah 



' ah 



M{r ah ) 



(III. 10) 



with v a h — v(r a h(t),t). Inversely, we can also write it as 
a curve t a h(r), which represents the time at which the shell 
labeled as r becomes trapped. 

In order to understand when the central singularity may be 
visible to faraway observers in the spacetime, we shall analyze 
the time curve of the apparent horizon which is given by, 



t a h(r) = t s {r) 



v ah {r) 



M 



be 2A - l 



-.dv , (III. 11) 



where t s (r) is the singularity time curve, whose initial point 
is to = t s (Q). Near r — 0, equation (IHI.llb can be written in 



the form, 



tah(r) = t ah (0) + 



dt a h 



dr 



r + o(r 2 ) 



r=0 

,2\ 



= t Q + X i(0)r + o(r 2 ) 



(III. 12) 



From this, it is easy to see how the presence of stresses affects 
the time of the formation of the apparent horizon. In fact, all 
the initial configurations that cause xi(0) t0 be positive force 
the boundary of the trapped surfaces to form after the occur- 
rence of the singularity, thus leaving the 'door open' for the 
null geodesies to possibly come out, away from the singular- 
ity. 

We can verify the visibility of the central singularity 
through the behaviour of radial null geodesies. This means 
that in order to verify the visibility of the singularity we must 
check if there are any future directed null or non-spacelike 
geodesies that terminate in the past at the singularity (for a 
detailed analysis of outgoing geodesies in singular spacetimes 
originating from collapse see 11811 ). It can be shown that if the 
singularity is radially censored then it must be completely cen- 
sored. On the other hand, if there are any future directed non- 
spacelike curves coming out from the singularity and reach- 
ing external observers then the singularity is naked. Further, 
in |Q]] it was proven that xi(0) > is a necessary and suf- 
ficient condition for outgoing radial null geodesies to escape 
the singularity, at least locally. Outgoing radial null geodesies 
are given by the equation, 



_ e i>-u _ R C ~ V ~ A 
dr ' y/b 



(III. 13) 



with the positive sign considered in order for t to be increasing 
moving away from the singularity. 

Now, introducing the variable ti = r™, we can write equa- 
tion dHI.l 31 ) as 



dR 1 



dv 



r a — l 



R' + Re^- V 



ar u 



VbVR 



'l_i4l/£ + (V I i .(III. 14) 



We thus look for a radial null geodesic terminating in the past 
at the singularity, whose equation near r = takes the form 



R = xqu, 



(III. 15) 



with xo finite and positive. Such null geodesies terminate at 
the singularity with a definite tangent and correspond to the 
curve 



t(r) =t + x a r a 



(III. 16) 



For a specific choice of a = |, we can solve equation ( 1III.14I) 
and evaluate the value of the tangent to the null geodesies at 
the central singularity which is given by r = 0, t = to- Thus 
obtaining, 



R 



xo = lim lim — 



dR 
du 



(0,to) 



^v/MoXi(0)V 



(III. 17) 
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which is positive and finite at the central singularity if xi (0) > 
0. It is now clear how xi is directly related to the tangent to 
the null geodesies at the singularity. In this case we conclude 
that the future directed radial null geodesies do come out of 
the singularity to reach far away observers. 

Let us stress that it is the sign of the quantity xi(0) that 
uniquely determines the visibility of the singularity, as in 
fact Xi(0) > implies that the time curve of the apparent 
horizon is future directed increasing, and future directed null 
geodesies in that case do come out. Therefore the central sin- 
gularity occurring at to is not covered. 

All those initial configurations for which M'(0), 6q(0) and 
a\(v) cause Xi(0) to be positive will force the apparent hori- 
zon to appear after the formation of the singularity, and will 
allow the null geodesies to come out of the singularity that 
forms at the center of the collapsing cloud. We conclude that 
the condition xi (0) > not only implies that the trapped sur- 
faces form at a later stage during collapse than the singularity, 
but it is also the sufficient condition for null geodesies from 
the singularity to be visible (at least locally) to external ob- 
servers. 



We see now immediately that the behaviour of the tangen- 
tial pressure (be it positive, negative, increasing or decreasing) 
is here reflected in the terms a(0, v) and a'(0, v) in xi (0), and 
can influence the final outcome of the gravitational collapse. 
In fact, it is the choice of the coefficients of the initial den- 
sity, velocity, and stress profiles that determines the quantities 
M, b and a, which are responsible for the behaviour of xi (0), 
which in turn determines the tangent to the singularity curve 
at the origin. As we noted earlier, the sign of xi(0) uniquely 
determines the final fate of collapse in terms of either a black 
hole or a naked singularity. As it was shown in fljj] , positivity 
of xi (0) is a sufficient condition for outgoing null geodesies 
to come out from the singularity, thus making it visible to ex- 
ternal observers. In some cases of physical interest it turns 
out that xi(0) = 0. This is likely to happen when marginally 
bound collapse with only quadratic terms in the expansion of 
energy density and pressure is considered, but it is not the only 
possibility. In general, whenever xi(0) = 0, the analysis of 
higher order terms must be carried out. In this case, we can 
write t s (r) as 



t s (r) =io+r 2 X2 (0) + o(r 3 ) 



(IV.6) 



IV. SMALL PRESSURE PERTURBATIONS IN DUST 
COLLAPSE 



and evaluate X2(0) 



d t 

dl'- 



as 



We would now like to examine here, using the formal- 
ism above, how the final fate of the gravitational collapse 
is affected when pressure perturbations are introduced in a 
pressure-free inhomogeneous dust cloud, which was other- 
wise going to terminate in a black hole final state. This pro- 
vides a useful insight into the role of pressure in gravitational 
collapse dynamics. Certain classes of collapse models with a 
non-zero tangential pressure are analyzed, their properties are 
discussed, and the collapse end-states are examined. 

For the sake of definiteness, we consider the case when the 
Misner-Sharp mass is given by an expression near the center 
as, 



M(r) = M + Mir + M 2 r 2 + o(r 2 ). 



(IV. 1) 



Further expressing the function g(r,v) (and therefore the 
tangential pressure) in the form, 

g(r,v) = g (v) + gi{v)r + g 2 (v)r 2 + o(r 2 ) , (IV.2) 

we have, 

r 2 ( 9 \ 

Pe = ( 3M .9o + 4Af lff0 r + -M o9l r + ... j . (IV.3) 

With these expressions, we are then able to write the quantity 
Xi(0) in dependence of the physically relevant profiles M(r), 
b(r), a(r, v) as, 

Xi(0) = -- / — w dv , (IV.4) 

2Jo (^+6 (0)+2a («)) 3 

where the function a(r, v) (which is related to v via equation 
(III. 14b ) has been written as, 



X2(0) 



(2«i+^(0)) 2 _ , 
(^+6 (0) + 2a ) i 

1 j& + 2a 2 + 2a 2 + 2& (0)a + ^(0) 

2 (Mo +6o (0) + 2a o ) 1 



9o(v) 



Mo 



6 (0) + 2a 



dv 



(IV.7) 



We also note here that the negative tangential pressures are 
not forbidden by the energy conditions. Here we observe, how 
the presence of negative pressures, that may be due to the pres- 
ence of exotic matter in the form of dark matter or some other 
repulsive effects, may influence the final stages of collapse, 
since it will affect the value of xi through the function a\(v) 
given in ( III. 221 ). 

The second order term in the equation for the singularity 
curve is given by equation ( IIV.7I ). It turns out that this term 
becomes relevant in many physically realistic scenarios which 
consider the mass profile and the density profile to be written 
with only quadratic terms in r. From the regularity conditions 
analyzed earlier, we can further see that M'(0) = (that is 
reflected in the requirement that M\ = 0) is a reasonable as- 
sumption. This is consistent with the analysis of energy den- 
sity profiles, that, for physically realistic scenarios, have only 
quadratic terms in the expansion. The energy density p, can 
thus be written as p(r,t) = pa(t) + p2(t)r 2 + p±(t)r A + .... 

In the situation where mass and density profiles are given 
only in quadratic terms, we then have, 



a(r, v) — a (v) + a\{v)r + a 2 {v)r 2 + o(r 2 ) . (IV.5) 



M 2n +1 — d2n+l — #2n+l — 



(IV.8) 
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We can now examine a few physically relevant classes of 
collapse models with such a property, together with the pre- 
scription for marginally bound collapse (60 = 0). This choice 
of the mass, pressure profile and velocity profile causes xi (0) 
to vanish and the final outcome of collapse is decided by the 
X2(0) term, which in this case reads, 



X2(0) 



1 + 2a 2 + 2a 2 + 2g Q + 2a ) 



-Mo 



2a o y 



dv. 



(IV.9) 

The analysis is then carried out along the same lines described 
above. 

Firstly, an important class of collapse models, also to be 
noted as the reference frame, is that of an inhomogeneous dust 
cloud. Inhomogeneous dust has been widely studied and the 
occurrence of naked singularities in such cases is well-known. 
From the framework given above, we recover the inhomoge- 
neous dust when a(r, v) = g(r, v) = 0, and that together with 
the requirement that only quadratic terms appear in the den- 
sity, yields 



dust 
A2 



(0) = 



M 2 Jv , 1 M 2 

~ dv — ~o T 

M 2 6 M 2 



(IV. 10) 



The collapse will end in a naked singularity for all nega- 
tive values of M 2 . Further, in this case, the structure of the 
apparent horizon can be evaluated explicitly. Note that for 
M(r) = Mo the above case reduces to the Oppenheimer- 
Snyder-Datt collapse scenario, describing homogeneous dust, 
and all the terms x«s vanish as the singularity is simultaneous. 

Let us now consider the effect of introducing small tangen- 
tial pressure perturbations in an otherwise pressure-free inho- 
mogeneous dust cloud collapse, that was going to a black hole 
final state. To this aim, we choose a class of collapse mod- 
els restricted by the assumption that the second order term 
in the mass profile vanishes. We see that if M 2 — then 
xi ust (0) = 0> an d the fi na l outcome of collapse for the inho- 
mogeneous dust cloud is then decided by the next order term, 
namely X4 US '(0) = — \-^jr- Exactly in the same manner as 

described above, all positive values of M4 will cause the col- 
lapse to end in a black hole. If now we introduce an arbitrarily 
small tangential pressure of the form (IIV.3b . we see that near 
r = this can cause the collapse to end in a naked singu- 
larity, whenever the function go(v) is chosen in such a way 
that X2(0), as given by equation JIV.9I > with M 2 = 0, is posi- 
tive. We therefore have constructed a class of models of small 
tangential pressure perturbations of the LTB collapse scenario 
that drastically changes the final outcome of collapse on intro- 
duction of small pressures. 

A second illustrative class of models is that of a 'quasi- 
Hookean' pressures, which is given by a specific choice of the 
tangential stresses corresponding to 



(-l)" + Vo 



1 



(IV. 11) 



for n = 1,2,3.... This choice leads to an energy function 



G(r, v) of the form 
G(r,v) = -— 



1 + b r 2 



Mo 



1 



l+r- 



V s ) 



(IV. 12) 



In this case, putting again b — 0, we get 



1 f 1 



M 2 v 2 + 2M /i 



(M v + 2[Iqv 2 — 2/i 
1 ypv (v 2 -l) (! + £)) 

(M v + 2fi v 2 -2// )* 



-d«+(TV.13) 



dv . 



Some more considerations on the dynamical behaviour of this 
particular case will be presented in section [VTl also. What we 
can see from the above is that the values of the constants M 2 
and /io are responsible for the final outcome of collapse in 
terms of black hole or naked singularity, and that either of 
the outcomes are possible depending on the choices of these 
parameters. 

Finally, it is interesting to study the behaviour of a subclass 
of the previous model, choosing the class of small tangential 
stresses in the form, 



R M g 2 4 6 
Pe = 6 — ^— r + o(r°) . 



(IV. 14) 



In this case we see that ao = a\ = 0, and near r = we have 
A(r,v) — a 2 (v)r 4 + o(r 6 ) and g 2 — ja 2 _ v . Then the final 
outcome of collapse is decided by the sign of 



xa(o) 



dust 
X2 



1 07(M 2 + 2a 2 v) 



dv 



(0) ~T f 1 

M n 2 JO 



a 2 v*dv . (IV.15) 



We see that in this case it is possible that the introduction of a 
suitable tangential stress in the inhomogeneous dust cloud will 
cause it to end in a naked singularity, while the same model 
with a vanishing pressure resulted in a black hole, and vicev- 
ersa. The condition to be satisfied by the tangential stress is 
3 J a 2 vzdv < —M 2 . 

To see this more explicitly, in analogy with the previous 
example, consider for instance 



a 2 = -Mo 



1 - 



(IV. 16) 



then the tangential stress near the center of the cloud will be 
given by, 



Pe 



3 — r +o(r D ). 



(IV. 17) 



We note that the tangential pressure as well as the density di- 
verge at the singularity. The collapse here ends in a naked 
singularity whenever 



M 2 < 



24 



-Mo 



(IV. 18) 
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FIG. 1 : The Naked singularity region and the black hole region in 
the parameter space {M2,pe} for fixed values of Mo, r and v. 



as can easily be calculated from equation (IIV. 15b - 

In Fig. [T| the occurrence of black hole and naked singular- 
ity phases in the latest example are shown, depending on the 
possible values chosen in the parameters space of pe and M2, 
for fixed values of Mo, r and v (taken at the initial time). The 
line X2 = corresponds to pg = — | M °^ /2 r 4 , and acts as 
a critical surface separating the two outcomes, in which case 
the evaluation of higher order Xi is necessary in order to deter- 
mine the final fate of collapse. The axis pg — corresponds to 
the LTB inhomogeneous dust case, with the origin being the 
Oppenheimer-Snyder-Datt model. 

We see clearly that to have a naked singularity in the dust 
case, we must require M% < 0. This is not the case anymore 
in the more general situation with non-vanishing tangential 
stresses. In fact, such scenarios allow for positive values of 
Mi to lead to a naked singularity, provided that the tangen- 
tial stress (and hence /zq in this example) is negative, as well 
as negative values of M2 leading to black holes whenever a 
suitable positive tangential pressure is introduced. The inter- 
esting feature here is that the introduction of tangential pres- 
sures does not rule out the possibility of formation of naked 
singularity. 

This shows that the conjecture suggested in ifioll does not 
hold (even in the more restrictive case of quadratic density 
and pressure profiles), and that naked singularities are a gen- 
eral feature of gravitational collapse with pressures. The con- 
jecture stated that given a known Lemaitre-Tolman-Bondi col- 
lapse model ending in a black hole, all quadratic pressure per- 
turbations to the model will necessarily originate a black hole. 
Now we see that adding a suitable generic pressure perturba- 
tion to the LTB model can in fact change the outcome of col- 
lapse. 

^From a physical point of view, the most relevant quadrant 
in the model shown in Fig. Q]is the first one (where pe > 



and M2 < 0), since it corresponds to positive pressures and 
radially decreasing energy density profile. We see that both 
the black hole and naked singularity outcomes are possible in 
this case. 

We shall stress however that the model illustrated above, 
while presenting a simple and intriguing structure, does not 
depict the only possible structure of collapse with tangential 
pressures. Other examples, with different choices of the func- 
tions M, b and a, might provide an entirely different picture 
for the final stages of collapse, thus indicating once more how 
rich and complex is the description of collapse resulting from 
general relativistic analysis. 

This model shows also how the black hole formation pro- 
cess described in the well known homogeneous dust scenario 
is not 'stable' when arbitrarily small stress perturbations are 
introduced. Since, as shown in Fig. [T] it lies on the critical 
surface separating the 'naked singularity configuration space' 
from the 'black hole configuration space', any close neigh- 
borhood of the OSD model will contain initial data leading to 
both outcomes. 

On the other hand, we see that in the context of this model, 
inhomogeneous dust collapse leading either to black hole or 
naked singularity (thus away from the critical surface where 
X2 = 0), is 'stable' under small pressure perturbations, in the 
sense that it is always possible to find a small neighborhood 
of the configuration space which will develop the same final 
outcome. 

The above considerations illustrate the role pressures can 
play towards determining the final fate of collapse. Clearly, 
depending on their nature, the non-zero tangential pressures 
can help create a naked singularity, or a black hole. It follows 
that a mere introduction of pressure in the collapsing cloud, 
by itself, cannot help us restore the cosmic censorship. Finally 
we note that the structure illustrated for the above model ex- 
tends naturally to more general pressure profiles, still leaving 
unchanged the general behaviour where both naked singulari- 
ties and black holes might generate at the end of collapse. 



V. GENERICITY OF BLACK HOLE AND NAKED 
SINGULARITY FINAL STATES 

The consideration above helps us gain an important insight 
into the nature and genericity of the black hole and naked sin- 
gularity formation in gravitational collapse. Let us take, for 
example, a scenario where the term \2 is positive, thus lead- 
ing the collapse to a naked singularity final fate. In that case, 
it is clear that there always exists a small enough perturbation 
in the initial data, which is in the form of small changes in 
the initial density, pressures, and velocities of the collapsing 
shells, which will still preserve by continuity the positive sign 
for the xi term, thus preserving the naked singularity final fate 
of the collapse. A similar conclusion will hold for the black 
hole final state of collapse as well. 

In this sense, if a collapse evolution was to develop in a 
naked singularity, there is a small neighborhood around that 
specific initial set in the configuration space of initial data, 
which also continues to take the collapse to a naked singular- 
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ity. In such a sense, the naked singularity formation (as well 
as the black hole formation) is generic in the gravitational col- 
lapse with non-zero tangential pressures case. 

In the light of these results, we can see that the idea that 
the introduction of pressures must favor the black hole forma- 
tion process is not correct, at least for non-vanishing tangen- 
tial stresses. One may argue that this situation arises because 
we have neglected radial pressures, a choice that, while still 
not being the dust model, is not sufficiently 'realistic' . Nev- 
ertheless the same arguments can be extended to the case of 
generic pressures (as it will be done in a future work), to show 
that collapse scenarios leading to a naked singularity have a 
small neighborhood in the configuration space of initial pres- 
sure profiles, that still leads to a naked singularity. 

Since the parameters that determine the final fate of col- 
lapse cannot be chosen arbitrarily if one wishes to obtain a 
black hole, we can say that a specific final outcome of col- 
lapse requires some tuning in the initial data in contrast with 
the possible view that any initial configuration with pressures 
will necessary destroy the naked singularity formation picture 
and lead to a black hole. 

The fact that clearly emerges from the present analysis is 
that the scenarios describing the formation of naked singular- 
ities as the endstate of collapse are not mere counterexamples 
to cosmic censorship whose relevance is limited and confined. 
We can therefore say that the presence of tangential stresses 
can 'undress' a singularity that was covered in the absence of 
stresses, or viceversa, a naked singularity dust model can be 
covered by the introduction of a suitable positive pressure. 

It is therefore useful to relate the scenario presented here 
with the known collapse models without pressures, namely 
the dust collapse scenarios. The first collapse model within 
the framework of general relativity was studied by Oppen- 
heimer, Snyder and Datt, and it described the collapse of a 
massive sphere made of homogeneous dust. The Lemaitre- 
Tolman-Bondi models then considered the collapse where in- 
homogeneous dust is allowed. It turned out that when inho- 
mogeneities in the density distributions are allowed, then the 
appearance of naked singularities takes place as possible end- 
state of gravitational collapse. Dust collapse models are well 
known and widely studied and the structure of the formation 
of trapped surfaces and of the singularity is well established, 
also in the case of higher spacetime dimensions lfl9ll . 

The dust collapse is obtained from the above framework, 
when we impose p r = pg = 0. In that case, from equation 
( III.5I I it follows that v' = 0, which together with the condition 
i/(0) = gives v — 0. This in turns implies that G = b(j-) 
and H = R 2 . For inhomogeneous dust in general we have, 



M'(0) 



-dv . 



(V.l) 



M(0) 



Mo) 



We see that inhomogeneous dust can lead to the formation of 
naked singularities or black holes, depending on the behaviour 
of the mass profile and the velocity profile 112011 . Furthermore, 
in the marginally bound case it is possible to evaluate explic- 



itly the apparent horizon from equation ( Mil 11 ) 



as 



1 Mi I M? 1 M 2 

tah{r) = t - - — T r + - — j- - - — 3 

6 M* V 4 Mn 6 M n 2 



r 2 + o(r 3 ) 



(V.2) 

while in the bound and unbound cases, the apparent horizon 
curve can be read from equation ( MI 111 ). 

Collapse is said to be homogeneous when p = p(t). Ho- 
mogeneity, together with the initial condition v{r,tj) = 1, 
implies some constraints on the possible matter models and 
velocity profiles. Now let us consider the following three con- 
ditions, 

(i) M = M a , 

(ii) v = v(t), 

(iii) bo(r) = k, 

where the second condition implies a simultaneous collapse, 
in which all shells of matter fall in the singularity at the same 
time. We can see that homogeneity, dust and the above con- 
ditions are not all mutually independent. In fact, in the case 
of collapse with vanishing radial stresses, which is under con- 
sideration here, we can prove the following statements: 

1 . Homogeneity (i) and (ii), 

2. (ii) <f^ (i), (iii) and dust. 

We can see from this that homogeneous collapse implies 
necessarily dust, while vice versa is not necessarily true (note, 
however, that when radial stresses are included it is possible to 
have homogeneous collapse also in the case of a perfect fluid 
with homogeneous pressures). Therefore, the Oppenheimer- 
Snyder-Datt model implies necessarily that the three condi- 
tions above are satisfied. We note that the simultaneous col- 
lapse (which implies homogeneous collapse in the case when 
M = Mo) implies that all shells fall into the singularity at the 
same time t = to, thus excluding the possibility of its visi- 
bility due to the structure of the apparent horizon in that case. 
The quantity Xi(0) vanishes identically in this case for every 
i. In order for the singularity arising in dust collapse scenar- 
ios to be visible, inhomogeneities in the density distribution 
must be allowed for and have to be taken into consideration, 
thus dropping the condition (i), which automatically implies 
the non-occurrence of a simultaneous collapse. 

To see the first statement, it is sufficient to use equation 
(111.41 ) with the requirement of homogeneity. From that it fol- 
lows that v(r,t) — 3CM(r)a(t), which together with the 
requirement that v(ti) = 1 implies (i) and (ii). To see the 
inverse implication, it is sufficient to put (i) and (ii) into equa- 
tion ( 111.41 ). To prove the second statement, we make use of 
equation (111.7b written in the form, 



M(r) b(r)e 2A - 1 



h(t) 



(V.3) 



Since h does not depend on r, it has the same value for any r, 
which is also that obtained for r = 0, namely, 



h = 



Mo 

v 



•Mo) 



(V.4) 
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^From this it follows immediately, 

M(r) = Mo, ff (r,«) = a(r,v) = 0, 6o(r) = 6 (0) . (V.5) 

Vice versa, assuming dust, and conditions (i) and (iii) to hold, 
we see that v does not depend on r, which implies v(r, t) = 
a(t) + (3(r). That together with the requirement v(ti) = 1 
implies v = v(t). 



VI. GRAVITATIONAL COLLAPSE WITH 
COSMOLOGICAL CONSTANT 

It is interesting and indeed important to see how the pres- 
ence of a non-zero cosmological constant can affect the dy- 
namical evolution of a gravitationally collapsing massive mat- 
ter cloud. The interest in models with cosmological constant 
arises primarily from the observation of the accelerating speed 
in the expansion of the universe, an observed phenomenon 
that indicates that the universe might be filled with some re- 
pulsive force known as 'dark energy'. One of the possible 
ways by which dark energy can be described is through the in- 
troduction of a non-zero cosmological term in Einstein equa- 
tions. 

The earlier works on the subject were constrained by some 
simplifying assumptions. Namely the collapsing matter was 
taken to be a dust cloud or the function v was chosen to have 
a specific dependence on t and R in order to simplify the in- 
tegration of equation ( III. 131 ) (see l2llD . 

In the present section we extend the previous results to the 
most general case and test how the dynamics and the final 
stages of collapse of the models described above are affected 
by the introduction of the cosmological constant. In doing so, 
we consider an example of dynamical collapse in the pres- 
ence of tangential stress and cosmological term, that connects 
directly with the models described in section HVl 

The reason for the interest in the introduction of only tan- 
gential stresses in the collapsing cloud with cosmological con- 
stant resides in the fact that, from equation ( III. 191 ). we see 
that the matching with an exterior vacuum background must 
be done with the well-known Schwarzschild-deSitter or anti- 
deSitter spacetimes. On the other hand, the presence of ra- 
dial pressures would require the matching with a more com- 
plicated generalized Vaidya metric. 

Einstein equations ( 111.31 ) and ( III. 41 ). in the presence of a cos- 
mological constant term take the form, 



A 



P + A 



F 

RiR 
F' 



(VI.l) 



(VI.2) 



R 2 R' ' 

while the other equations remain unchanged. The first one 
above leads to a Misner-Sharp mass of the form 

1 



F(r,t) = r 3 M(r) + -AR 3 , 
3 



(VI.3) 



while the energy density remains unchanged. The energy con- 
ditions here imply 

P>0, P + A>0, p + pg >0, (VI.4) 



however, all of the above need not be respected now, and 
we may have a weaker form of energy conditions holding, 
depending on what the sign of the cosmological constant is 
chosen to be. The matching with the exterior metric in this 
case needs to be done with a spacetime that is asymptotically 
either Schwarzschild-deSitter or Schwarzschild-anti-deSitter, 
depending on the sign of the A term. 

The integration of Einstein equations proceeds exactly in 
the same way as described above, and we therefore obtain, 



v 2 



-v i'M 1 . a 
' h -Av 2 



be 



2A 



1 



= V(r,v), (VI.5) 



where the quantity V(r, v) can now be interpreted as an ef- 
fective potential. The regions of allowed motion are those for 
which V(r, v) > (since v 2 > 0). The zeros of V corre- 
spond to turning points of the dynamics where the collapse is 
halted and the shell bounces back. Therefore, if V(f, v) = 
for some f, the cloud will bounce back and not collapse in- 
definitely towards the center. Of course, this situation is not 
allowed in the dust cases when A = 0, r is small and where 
gravity is the only force acting on the collapsing particles. 
Still, bouncing behaviours were found to be possible in sce- 
narios with tangential stresses and bound velocity profiles. In 
the same manner, whenever A ^ 0, it is possible to observe 
such bouncing behaviours in the dynamics of the collapsing 
shells. 

Typically, such a scenario will become relevant for collaps- 
ing clouds of rather large sizes in the universe, for example, 
for the case of a large cluster of galaxies collapsing under its 
own gravitational attraction, or possibly for a large collaps- 
ing star cluster. In any case, for the sizes of a collapsing star 
of an average size, the effects of a non-zero A may not make 
any significant impact on the final outcome of its gravitational 
collapse. 

Following a similar analysis of the time curves as that de- 
veloped in section[ni] we can evaluate t(r, v) in this case as, 



t(r,v) 



zdv 



— + iA«2 

v 3 



be 2A -l 



(VI.6) 



This leads the numerator of \i (0) (which is the essential el- 
ement in determining its sign) to have exactly the same value 
as in equation dIV.4l >: 



Xi (0) 



M j L+b> (0) + 2ai(v) 



-dv 



(Mo + 4 w a + 6o (o) + 2ao(t;)) a 

(VI.7) 

We can thus see that the presence of the A term affects col- 
lapse in the sense that it might introduce some positive roots in 
the effective potentials, and therefore would give rise to some 
bouncing behaviours. But it will not affect the terms that are 
responsible for the sign of Xi(0)> which determines the final 
outcome of collapse to be a black hole or a naked singular- 
ity. This appears reasonable, when we think of the presence 
of the cosmological constant as related to a dark energy in 
the universe, because the effect that the dark energy will have 
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on collapse can be that of disrupting and preventing collapse 
from happening by causing a bounce. But it cannot influence 
the nature of the final fate of a collapse that is uniquely deter- 
mined by the initial configurations for densities, pressures and 
velocities. 

Furthermore, the presence of an added positive A-term can 
lead to M + \\ v 2 + > while — + *s^=i < 0, 

v 3 v r z ' 

thus in fact allowing the evolution for some configurations that 
were forbidden in the case of vanishing A. 

The visibility of the singularity in the case of collapse with 
a cosmological constant differs from the case studied above 
also in the structure and formation of the apparent horizon. 
From equation ( IVI.3I I and equation ( IIII.9I I we get 



We therefore choose: 



' ah 



M 

Vah 



A 



' ul, 



(VI.8) 



from which we see that in general the apparent horizon equa- 
tion is a cubic equation in v that can have zero, one or more 
than one positive roots, thus affecting considerably the forma- 
tion of trapped surfaces. Still it is possible to show (see I2H0 . 
that the presence of the cosmological constant does not pre- 
vent the possibility of formation of naked singularities, which 
continues to be a general feature of gravitational collapse. 

It is interesting at this point to analyze some examples of 
dynamical behaviours based on the above considerations. We 
will thus make a specific, although reasonable, choice of the 
mass, pressure and velocity profiles and confront the effective 
potential given by equation ( IVI.5b written in the form, 



V(r,v) = 



M W , 1 . A « 2 + - (r ' u)_1 



(VI.9) 

for different values of A. 

Following the examples provided in section [IV] we will 
choose only quadratic terms in the mass and stress profiles 
and study a 'quasi-Hookean' equation for the pressure of the 
form 



-A(r,v) 



= 1 - HO 



1 + r 2 



1 - 



1 



(VI. 10) 



It is easy to check that in this case ao = and a,2 n -2 = 
E^UC - l) n r 2n (n = 2,3...), and we obtain a model that 
has the same singular behaviour of that analyzed in equation 
( HV. 15b . Therefore, for a suitable choice of the parameters, ac- 
cording to the inequality given by ( II V. 18b . the final outcome 
of such a scenario can be either a black hole or a naked sin- 
gularity. Furthermore, in order to keep the example as gen- 
eral as possible while maintaining the mathematics behind it 
manageable, we will consider bound, unbound and marginally 
bound collapse models only in the case of 60 constant. Finally, 
for simplicity, in the mass profile we will consider only one 
non-vanishing term in the expansion. This choice does not af- 
fect the dynamical behaviour very much, since the lower order 
terms in the expansion are the ones responsible for the general 
features of the mass profile in a neighborhood of the center. 



M(r) = M + M 2 r 2 , 
b(r) = 1 + kr 2 , 

1 + kr 2 

G(r,v) ■■ 



1 - Hoj^pi (1- ^r) 



(VI. 11) 
(VI. 12) 

(VI.13) 



As seen before the following model can lead to a black hole 
or a naked singularity, depending on the values chosen for 
M2, Ho an d k. Some plots of the dynamics of this model for 
a specific choice of the parameters are given in the figures. 
The functions so defined correspond partially to the functions 
analyzed by Magli in [10] (where the cosmological term was 
not considered) and appear reasonable in the sense that the 
pressure is small compared to the energy density. 




FIG. 2: Sample plot of V(r, v) in the unbound case, at a fixed value 
of r close to the center, with M = 0.01, M 2 = -0.001, ju = 0.1 
and k = 1. 

The first thing that can be noted is that the introduction of 
a positive cosmological constant can lead to unbound mod- 
els even in cases that were bounded in the absence of the A- 
term, thus making the dynamics possible for configurations 
that were previously not allowed. On the contrary, a negative 
cosmological constant term might introduce zeroes in V(r, v), 
thus altering and limiting the range of allowed configurations. 

We can see that, since the dynamics happens in the region 
of positive V, when V has no zeroes then no bouncing be- 
haviour can happen. Since in the region close to the center V 
diverges to plus infinity, we can infer that the core of the cloud 
will always collapse to the singularity. On the other hand, we 
note that at larger values of r oscillating behaviours between 
two fixed radii are possible. Furthermore at larger values of 
r bouncing behaviours can happen whenever the shell labeled 
with r reaches an event v for which V(r, v) = 0. In this case 
the collapse halts and the dynamics is reversed. In the figures 
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FIG. 3: Sample plot of V(r, v) in the marginally bound case, at a 
fixed value of r close to the center, with M = 0.01, M 2 = -0.001, 
/io = 0.1 and k = 0. 




FIG. 4: Sample plot of V(r, v) in the bound case, at a fixed value of 
r close to the center, with M = 0.01, M 2 = -0.001, fi =0.1 and 
k = -0.1. 



we have restricted the analysis to a close neighborhood of the 
center of the cloud. 

As we have noted earlier, the key physical interest for the 
bouncing behaviour comes for the collapsing clouds which 
are large enough in physical scales, but for stellar collapse, 
the cosmological constant value plays no significant role. Of 
course, even without a cosmological term, a non-zero pressure 
can cause a bounce in the outer regions of the cloud. However, 
the inner region must still collapse to the singularity necessar- 



ily, where the behaviour of the cloud is very much like the 
dust collapse. This is consistent with the regularity condi- 
tion we mentioned earlier, namely that the tangential pressure 
must go to a vanishing value in the limit of approach to the 
center. A negative sign for the cosmological term preserves 
the same qualitative nature of the collapse as above, however, 
a positive cosmological constant brings in interesting changes 
in the collapse dynamics, as noted above. 



VII. CONCLUDING REMARKS 

The framework to describe collapse of massive objects in 
general relativity is extremely rich and complex. Even in the 
limited case of spherical symmetry a wide array of different 
situations and collapse scenarios arises. 

We have here presented a simple and straightforward anal- 
ysis of the gravitational collapse of a spherical distribution of 
mass in the presence of tangential stresses. We have shown, 
how the occurrence of trapped surfaces and the eventual vis- 
ibility of the singularity forming at the end of the collapse 
process is directly related to the nature and behaviour of the 
singularity curve. 

This result further supports the conclusion that for a generic 
spherical object, the end-state of complete gravitational col- 
lapse will be either a black hole or a naked singularity. The 
initial configuration as well as the equations describing the 
properties of the matter cloud during collapse are the only fea- 
tures that, once coupled to Einstein equations for the dynamics 
of the model, will determine the final fate of collapse. 

Since the radial and tangential pressures play different roles 
during the evolution of collapse dynamics, the present work 
helps understand how much the outcome is due to either of 
them. By investigating the role played by tangential pressures 
near the formation of the singularity in full generality, with- 
out any restrictions on the type of matter responsible for the 
pressure itself, we allowed for the possibility to consider nega- 
tive stresses also. These, while still satisfying the weak energy 
condition (i.e. not considering phantom sources, and therefore 
constraining the pressure to a lower finite bound), are possible 
candidates for those exotic matter sources that have recently 
become of great interest in both cosmology and astrophysics. 

In the past years many studies have been devoted to under- 
stand the physical processes involved in the complete gravi- 
tational collapse of a massive star. The treatment presented 
above offers a comprehensive tool to deal with the occurrence 
of naked singularities in the case of vanishing radial stresses 
and consequently provides a general picture in which earlier 
results are included. As an example, the work by Magli on 
collapse with tangential stresses [10], that provides a complete 
characterization of the outcome of collapse in terms of the root 
equation developed by Joshi and Dwivedi, can be translated in 
the present framework. As noted, the use of mass-area coor- 
dinates (or the area-radius coordinates) ultimately reduces to 
the change of variables (r, t) — > (r, v), made here in order 
to study the equation of motion (III. 18b . Also, the analysis of 
the Joshi-dwivedi root equation, which governs the visibility 
of the singularity, is equivalent to the analysis of radial null 
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geodesies exposed in section Hill 

Furthermore, the results that give the conditions regarding 
the strength of the curvature singularity occurring at the center 
which were obtained by Nakao et al. IU2I1 can be translated in 
the above formalism. Equations (2.13) and (2.14) in that paper 
are then found to be equivalent to equations (IIII.4b and (III. 18b 
which determine the equation of motion and the expansion 
of t(r, v) near the center to first order in r, while equation 
(2.18) gives a condition for the occurrence of outgoing null 
geodesies equivalent to equation (IIII. 1 7b . 

The models provided in section [IV] and section [VT] can be 
connected to the ones discussed in the previous papers and 
therefore provide a useful tool to broaden the spectrum of col- 
lapse scenarios studied in detail and expand our comprehen- 
sion of the final stages of complete gravitational collapse in 
general relativity, from the analytical viewpoint. 

The fact that positivity of the first non-vanishing Xi (0) is a 
necessary and sufficient condition for the occurrence of naked 
singularities may be possibly used in computer simulations, 
to analyze the endstate of collapse in a broader and more real- 
istic range of situations. In fact, the early works in numerical 
relativity by Shapiro and Teukolsky 02211 on naked singular- 
ity formation were hindered by the fact that the connection 
between the delay of trapped surfaces formation and the ex- 
istence of outgoing null geodesies was still unclear. Shapiro 



and Teukolsky considered models describing the collapse of 
prolate spheroids (thus in connection with the problem of axi- 
ally symmetric collapse and the hoop conjecture), and models 
describing the collapse of counter rotating particles with van- 
ishing angular momentum (the 'Einstein cluster' cited above). 
In both cases they found evidence for the non-occurrence of 
the apparent horizon before the time of the formation of the 
singularity. Still, Wald and Iyer showed that the fact that on 
certain slices trapped surfaces form at a later stage does not 
necessarily imply the visibility of the central singularity (they 
pointed out a time slicing of the Schwarzschild space-time 
whose evolution approaches arbitrarily the singularity with- 
out the appearance of an apparent horizon in any slice) II23I1 . 

We can argue here that in comoving coordinates the occur- 
rence of trapped surfaces and the visibility of the singularity 
are related via the sign of Xi(0)> thus providing a firm foot- 
ing for further numerical investigations. Hopefully, numerical 
simulations of naked singularity formation will help shed a 
better light on the final stages of realistic collapse and open 
new possibilities for the investigations of what happens in last 
moments of the life of a star. From an astrophysical point of 
view, this approach is clearly important, because if such pro- 
cesses indeed happen in the final stages of collapse of massive 
objects, it is then crucial to understand what kind of signature 
they bear in order to be able to eventually 'observe' them. 
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